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The modification of tlie Vlasov equation, in its standard form describing a diarged particle distri- 
bution in the six-dimensional phase space, is derived explicitly within a formal Hamiltonian approach 
for arbitrarily curved spacetime. The equation accounts simultaneously for the Lorentz force and 
the effects of general relativity, with the latter appearing as the gravity force and an additional 
force due to the extrinsic curvature of spatial hypersurfaces. For an arbitrary spatial metric, the 
equations of collisionless hydrodynamics are also obtained in the usual three- vector form. 
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I. INTRODUCTION 

The Vlasov theory for curved spacetime is well de- 
veloped, yet has been receiving but scant attention in 
plasma physics. The reason for this may be that the ex- 
isting theory relies on covariant formulation P, [2| , which 
renders difficult to use the intuition available through 
standard formulations of plasma physics. The so-called 
3 -|- 1 formalism could help solve this problem, by project- 
ing the general relativistic (GR) equations on the more 
conventional three-dimensional (3D) space; specifically, 
the Maxwell's equations and the particle motion equa- 
tions can be put in a three-vector form similar to that 
in the Minkowski metric However, a question re- 

mains how exactly the Vlasov equation, in its standard 
form describing a charged particle distribution in the 6D 
phase space is modified when spacetime is curved. 

To answer this question is the main purpose of the 
present paper. The same problem was previously ad- 
dressed ad hoc for a specific metric in Refs. [l^, [HI and 
more generally in Refs. [3 [T^. However, in the lat- 
ter case, electromagnetic interactions were not included, 
and the spacetime basis was not identified that would 
yield conventional three- vector equations of collisionless 
plasma hydrodynamics (see below). To close the theory, 
one thereby needs to specify the 3-1-1 equations of charged 
particle motion explicitly and in the form analogous to 
that commonly used for the Minkowski metric 

The routine approach to these equations that involves 
the full-fiedged machinery of differential geometry [6l-[l0| 
cannot be employed for this purpose. Also, their general 
form is reported by different authors in different forms 
which are neither manifestly equivalent, nor always ac- 
curate [l6[. Thus, it seems warranted to rederive the 
corresponding equations from scratch, particularly, in the 
Hamiltonian representation that is immediately applica- 
ble to the Vlasov theory. A spin-off here is that, once 
the Hamiltonian formalism is developed, the covariant 
approach used in Refs. [T^.IT5| becomes unnecessary. In- 
stead of generalizing Refs. [l4 [lH then, one can as well 
make the 3-1-1 Vlasov theory self-contained and refor- 
mulate it in a manner familiar from the standard plasma 
physics. This constitutes the second purpose of our pa- 
per. 



Finally, our third purpose is to derive hydrodynamic 
equations from the Vlasov theory in curved spacetime in 
their usual three-vector form. Unlike app roaches postu- 
lating a hydrodynamic closure [l^ [T7l - ll9{ . this will yield 
a fundamental fiuid treatment of plasmas in the collision- 
less limit, where the commonly used ideal-fiuid approxi- 
mation B i, [11 HI does not apply. 

Specifically, our results can be summarized as follows. 
The three-dimensional dynamics of a charged particle in 
an arbitrary spacetime metric, traditionally addressed 
within differential geometry, is reformulated in terms of 
linear algebra and Hamiltonian formalism. The modifica- 
tion of the Vlasov equation, in its standard form describ- 
ing a charged particle distribution in the 6D phase space, 
is then derived explicitly. The equation accounts simul- 
taneously for the Lorentz force and the effects of general 
relativity, with the latter appearing as the gravity force 
and an additional force due to the extrinsic curvature 
of spatial hypersurfaces. For an arbitrary spatial met- 
ric, the equations of collisionless hydrodynamics are also 
obtained in the usual three-vector form. 

The paper is organized as follows. In Sec.|ITl we restate 
the 3-1-1 formalism by amending the approach that was 
adopted in Ref. f^]. In Sec. Illli we derive the equations 
of individual charged particles interacting with an elec- 
tromagnetic field in space with an arbitrary metric. In 
Sec. IIVI we obtain the Vlasov equation in a number of 
equivalent forms. In Sec. |Vl we derive the equations of 
collisionless hydrodynamics. In Sec. IVI[ we summarize 
our main results. Supplementary calculations showing 
how our formalism relates to that in Ref. [1] are given in 
Appendix. 



II. SPACETIME GEOMETRY 

In this section, we restate the 3 + 1 formalism, which 
is based on the so-called slicing approach summarized 
in Ref. [3,] (not to be confused with the 1-1-3 formal- 
isrn, which is based on the so-called threading approach 
[S [HI)- Our purpose is to restate the known theory in 
a systematic, self-contained, and yet concise and simple 
form, relying on the reader's background in linear alge- 
bra rather than differential geometry; hence, the rules of 
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index manipulation, still required, are included here. For 
reviews and other formulations, see Refs. [l,iM[IMl. 

A. Metric tensor and index manipulation 

First, let us provide a brief introduction to the part of 
tensor analysis that will allow a reader, assumed familiar 
with linear algebra, to understand the rest of the paper 
(excluding the Appendix, which contains no new results). 
Suppose that the spacetime is equipped with a metric 
tensor g (the caret is to denote rank-two tensors) , which 
then defines a scalar product for any given four-vectors 
X and Y, 

X.Y^g(X,Y), (1) 

as a symmetric bilinear form. Consider a set of basis 
vectors e^, where the Greek indexes span from to 3. 
Then, each vector X can be decomposed as X = e^X'^ 
(and similarly for Y; summation over repeated indexes is 
assumed), yielding 

X.Y = 5^,X^y^ 5^, = g,^ = g(e,,e^). (2) 

Unlike in the Euclidean (or Minkowski) space, the met- 
ric coefficients 5^1, = • may not form a diagonal 
matrix, i.e., are not orthogonal to each other. Hence, 
additional, "dual" basis vectors are introduced [26i] via 

e'^ • e. = Sii, (3) 

which can also be understood as the definition of the 
"mixed" metric coefficients = g^'^ = S^. Then, the 
vector components X^^ in the original basis, or the so- 
called contravariant components, can be found as X^^ = 
• X. On the other hand, one may as well define the so- 
called covariant components as the vector components 
Xfj^ in the dual basis, as = • X; then, 

X^^e^-e.X" ^g^.Xr (4) 

Finally, let us define the inverse metric tensor g~^ with 
the metric coefficients denoted as g^" = g'^^^ (the symme- 
try being inherited from that of g^,y); hence, from Eq. (jlj, 

X" = g'^'X^. (5) 

Then, Eq. ^ can be equivalently put as 

X • Y = g^.X^Y'' = g'^'^X^Y, = X^Y^ = X^Y^. (6) 

In addition to the four-vector space, we will also be 
dealing with a three- vector space (Sec. Ill C|) . with objects 
to be denoted with bold (X instead of X). For those, 
the same rules of index manipulation apply, except with 
the metric tensor 77^- (instead of g^iv) and Latin indexes 
(instead of Greek indexes), spanning from 1 to 3. For 
further reading on tensor analysis in application to the 
general relativity, one is referred, e.g., to Refs. [27l - l29| . 



B. Spacetime basis 

Suppose that time t = (assuming the speed of light 
is equal to one) is defined as some function of space- 
time location, such that constant-t hypersurfaces E( are 
space- like. Introduce three arbitrary basis vectors ei{x'^) 
as tangent to these surfaces; hence the generalized coor- 
dinates in St, which we denote as space. Then, an 
arbitrary infinitesimal four- vector dx is decomposed as 

dx — Si dx^ + 60 di, (7) 

where eo is the basis vector along the time axis, yet to be 
defined. Since the four-gradient Vi is normal to (and 
thus orthogonal to e^), Eq. ^ yields 

Vt-dx^ (V< • eo) dt. (8) 

On the other hand, 

dt 

dt = - — dx^' = Vt ■ dx. (9) 

Then, from Eq. ([5]), one obtains 

Vt-eo = l, (10) 

meaning that is dual to eo, i.e., = e°. 
Consider a normalized vector 

n = -aVi, (11) 
where a is a scalar function such that 

n-n = -l, (12) 
and the sign is chosen assuming the metric signature 

(-,+,+,+) (13) 

(That is, n is the time- like unit normal to Ej.) Then, 
Eq. (fTU]) finally rewrites as 

n • eg = —a. (14) 

Equation (jl4p is the only requirement on how eo should 
be defined. Although the time axis could be normal to Ej 
(like n and Vt) , in general, eo can also have a component 
P tangent to Et. Thus, from Eqs. and ([11]), the 
general form of this basis vector is [Fig. [Ija)] 

eo = an + |3, (15) 

and the latter can be understood in two different ways. 
If eo is imposed, then functions a and |3 are used to 
'parameterize the given basis Bq = (eo, e^). Alternatively, 
one may be allowed to choose eo as needed; in that case, 
a and (3 are free parameters, and it may be convenient to 
pick them differently depending on a problem of interest. 
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FIG. 1: Schematic of the spacetime bases (not in scale), (a) Here n is the unit normal to the space hypersurface Et; eo is 
the basis vector that determines the time axis (t-axis); ei is the spatial basis vector (the other two spatial dimensions are not 
shown); C* = an determines ("-axis; a is the lapse function; (3 is the shift vector, (b) The vectors (eo,ei) form the basis So; also 
shown is the dual basis (6*^,6'). (c) The vectors (C*,ei) form the basis S^; also shown is the dual basis (C*,e'). 



To understand the physical meaning of the functions 
a and (3, consider the coordinate form of n in Sq: 



n" = e" • n = Vt • n = -(n • n)/a = (16) 

EE • n = [e* • (eo - (3)]/a = -^'/a, (17) 

no = eo • n = (on + |3) • n = —a, (18) 

rii = ei ■ n ~ —a (e^ • e'') = 0, (19) 

or, in a compressed form, 

= (1/a, -/3Va), = 0,) (20) 



[cf. Fig. IHa) and (b)]. Then, due to Eq. ([11]), one can 
treat n(a;'') as the four- velocity of some observer at x^, 
which we call, after Ref. [1], a fiducial observer (FO). 
Introducing the FO proper time Tpo, one gets [Eq. d^ ] 

a = dr^o/dt. (21) 

Therefore, a is called the "lapse function" . Similarly, /3* 
can be understood as minus the spatial velocity of FO, 
meaning that determine the rate at which the coordi- 
nate mesh on Et is shifting with respect to FO; thus, |3 is 
called the "shift vector" . Since the latter is not a physi- 
cal velocity, it can be arbitrary, including superluminal, 
and therefore eo is not necessarily a time-like vector. 

C. Spatial metric 

To switch from a spacetime, where an arbitrary vector 
X is decomposed as [see Eq. 

X = BiX^ + {an + ^) (22) 

to the three-vector representation, construct the spatial 
vector space as follows. Introduce a symmetric tensor 

h^^ + 71^^71,, (23) 



or h = I + nn (where I is the unit tensor), to project X on 
the plane tangent to Et: 

X=h-X = X+(n-X)n. (24) 

Since h • n = 0, directly from Eq. ((22)) one can see that 
X'' = {0,X' + P'X^). (25) 

For X is a tensor contraction, transform as vector 
components by definition, and thus so do X^ (but not 
necessarily X*). On the other hand, X'^ = 0; hence, 

X={X\X'^,X^) (26) 

can be considered as a spatial three- vector. 
Consider a length element in Sj: 

dx • dx = dx^dxi — dx^dXfj_, (27) 

where we used that dxp = 0. Notice further that 

dx h • dx, (28) 

and thus dx^ — /i^iydi". Then, Eq. ((27)) yields 

dx • dx = hf^iydx'^dx^ — rjijdx^dx-' , (29) 

where we introduced 77.^ = hij to distinguish the four- 
tensor hf^^ from its spatial part, which represents a sym- 
metric three-tensor. From Eq. (|29p . it is then convenient 
to address rjij as the spatial metric. One hence defines 
the mixed and contravariant metric tensors through 

^ = v'^Vk, = S). (30) 

This allows raising and lowering indexes as 

r^'^X, = X\ i^,,X' = X„ (31) 

where the covariant components satisfy 

X, = X„ (32) 

as flows from Eq. ((24)) and rii = 0. Hence, a three-vector 
scalar product can be defined as 

X • Y = Vv^'Y^ = ri'^X.Yj = X% = X,Y\ (33) 
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D. Spacetime metric 



The three- vector components, Eq. ((25|) . can be under- 
stood as components of X in the basis = (C*, e,), with 
the dual basis being (C*,e*), where 



an, 



(34) 



and e* = h-e' [Fig.[TIc)]. Denoting the vector components 
in B(; with underbars, one obtains, in particular, that 



(a \0'), = (-a,Oj), 



(35) 



and, for vectors X projected on the plane tangent to St 

[Eq. mi 



(36) 



In other words, the components of spatial four-vectors 
are the same in Bq and , and will not be distinguished 
from now on. Finally, from dx = an dt + dx [see Eq. ((22)) ]. 
we get 



dx ^ C* dt + dx, 



(37) 



so d( = C*'dx equals dt. (Yet, when taking partial deriva- 
tives, one must distinguish ^-axis, along which a;* are 
fixed, from x^'-axis, or i-axis, along which are fixed.) 
Hence, an arbitrary spacetime interval can be put as 



dx - dx ~ —a^dt^ + dx • dx. 



(38) 



Therefore, in the basis Bq, the metric tensor is repre- 
sented by a block-diagonal matrix 



-a2 
rjij 



(39) 



III. SINGLE PARTICLE MOTION 

In this section, we derive the motion equations for indi- 
vidual charged particles interacting with an electromag- 
netic field in space with an arbitrary metric "qij . We con- 
sider both the coordinate and the vector form of these 
equations. Particularly, this explains how Eq. (3.8) in 
Ref. [3], similar to those, can be obtained without using 
the concept of a covariant derivative (also see Appendix). 
Additional information on the Hamiltonian formalism for 
the particle motion in curved spacetime can be found in 
Ref. [si. 



A. Canonical equations 

Consider a particle with mass m and charge q inter- 
acting with an electromagnetic four-potential A, so the 
particle action S reads as (see, e.g., Refs. [30l Isij) 



S 



; \/—dx ■ dx + qA ■ dx 



(40) 



Introduce the particle Lagrangian L through S = J Ldt. 
Hence, from Eq. ([55]) . one gets 

L = -mx-^ + qAjiv' -f3') + qAn, (41) 

where the velocity components read as [Eq. (|25p ] 

= dxydt = v' + P\ (42) 

and, with v . v = , x is given by 

>f-i = Va2 - v2. (43) 

(As before, the bold symbol v denotes a three-vector 
with components u*, not = dx^ /dt.) Then, one puts 
Eq. (pij) in the following form: 



L = -m\/d^ - v2 + gv . A + qA^, (44) 
where Aj = Aj [Eq. and 

= C* • A - Ao ~ P'A.j (45) 

is the covariant component of A along ^-axis. 

The three canonical momenta P^, completing the 
canonical pairs (x^P^), are defined as Pi = dL/dv\ or 

Pi = mxT]ijV^ + qA,. (46) 

Using Eqs. (pij) and ([5^ . one obtains then: 

Pi — mxvi -f qAi — mxvi + qAi. (47) 

Thus, the three Pi given by Eq. ((46|) are equal to the 
covariant components of the corresponding four- vector P. 
Similarly, the kinetic momenta. 



equal 



Pi^ Pi - qAi 



Pi — mxrjijV^ — mxVi, 



(48) 



(49) 



and thus coincide with the corresponding components of 
the particle four-momentum p. Therefore, the formally 
introduced variables Pi and pi satisfy (unlike Pi and pi) 
the index manipulation rules introduced in Sec. Ill CI 
Hence, invert Eq. (|49p as 



(50) 



This allows the particle Hamiltonian, 

H = Piv' - L = o?mx - qAc^, (51) 
to be expressed as a function of pi. Specifically, 

>< = 7/a, 7 = yrTpVm^, (52) 
thus, iJ rewrites as 

H = amj — qA(^, (53) 
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where, in terms of the canonical momentum P, one has 

7 = ^1 + {V -qAf/m?. (54) 

Finally note that, since one can also put the latter as 

1 



where V is the particle velocity as seen by FO, 



V ; 



1 dx 



dt cIty 



(55) 



(56) 



one can understand 7 as the particle Lorentz factor as 
measured by FO. 

The canonical equations are now obtained as follows. 
For the coordinates, dx^/dt = dH/dPi yield the already 
known equation 



dx^ 
~dt 



ap 



(57) 



For the canonical momenta, one has dPi/dt = —dH/dx'^, 
or 



a dt ^ * 2mj a \ dx^ 



dx^ 



where we substituted 



dx^ 



(58) 



(59) 



[from differentiating Eq. (l30l) ] and introduced gj (not to 
be confused with the metric tensor g) as components of 
the three- vector 



g 



-Vln( 



(Unlike V, which is the gradient in the 4D space, the 
bold symbol V denotes the gradient in the 3D space.) 
Using Eq. (j48|) . one then gets for the kinetic momenta: 



l_ dpi_ 
a dt 



m7g 



dx^ 2m7 
with Ai being the Lorentz force: 



■A,: 



C 



dx^ 



dAj dA, 
dx^ dxi 



(61) 



(62) 



B. Lorentz force 

Consider expressing A^ in terms of 
_ dA, dA^ 



(63) 



which falls under the definition of a tensor (see, e.g.. 
Sec. 83 in Ref. [H or Sec. 4.2 in Ref. [H). By anal- 
ogy with the Minkowski spacetime, write this so-called 
electromagnetic tensor (Sec. 90 in Ref. [2^) as 

Fy,^ = n^,E^ - n^E^ 4- ext^v^B^n^. (64) 

Here E and B are four- vectors with zero time components 
(hence the bars), t\^uK is the permutation pscudotensor: 

<^\^lVK = \/^[A^j^k], (65) 

[XjjLVK] is the permutation symbol, and 

del g^u — —o?Tl, rj — detrjij, (66) 
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the equality flowing from Eq. (15^ . [This form ensures 
that Ffj_i, is indeed a tensor and introduces the exact 
amount of free parameters (the six nonzero components 
of E and B) to define an antisymmetric matrix like i^^j/.] 
In terms of F^j^, Eq. (|62p rewrites as 

A^ ^ ^ {F^c + K,v^) , (67) 

where f^i, are the corresponding components of i^^i, in 
the basis Bq. From Eq. (|64|) . one gets 

F^c^-ncE, + ex^aB'n\ (68) 

where we utilized Eq. (|36p for the spatial four-vectors. 
Employing Eq. p5p in the form 

(69) 



aEi 



Ui =0, Uc^ — —a, — a 
one gets 

FiQ = aEi + a~^€Qic^j& 
because [C^O] = 0- Similarly, 



Ftj = ex,jkB''n^ = a ^ 



(70) 



(71) 



Here we use e^fe = y^[ijfc] to define the permutation 
(60) three-pseudotcnsor eij^. Unlike Eq. (16) in Ref. [l^l: this 
is a standard definition of tijk^ which automatically en- 
sures that the Lorentz force has a three-vector form. 
Finally, Eq. ((67|) rewrites as 



K^^q{E^ + a h^jkV^B'' 



(72) 

With the definition ((56)) . one hence obtains that A^ can 
be regarded as covariant components of the three- vector 

A = g(E + VxB). (73) 



Equation ([73]) is similar to that in the Minkowski space. 
However, notice the difference between V that enters 
here and the velocity v = dx/dt [Eq. (157 ])]: in particular, 
see Appendix for comparison with Ref. [S^]. Notice also 
that, within the three- vector formalism that we adopt, 
the shift vector (3 does not explicitly enter the above 
derivation of the Lorentz force, unlike in Ref. 22|. It can 
be reintroduced, though, by substituting = + (5^ 
[Eq. (|42|)]. Contrary to Ref. H^, the so-called transport 
velocity = dx^ /dt here is also a three- vector, e.g., be- 
cause it equals the difference of three- vectors and . 
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Metric-caused forces 



Now let us revert to Eq. ([5T|) and calculate the effect 
due to the metric rjij. First, raise the index using 



P.- (74) 



Substituting Eq. ((57)) and emplo ying Eq. ([30)) . one ob- 
tains (similarly to Sec. 87 in Ref. |29l|): 



a dt 



TW-f 



■ TO7g 



a dC 



f+A', (75) 



where F^^ — F^^, are given by 



jk 



-11 



dxi dx'^ dx"^ 



(76) 



also known as connection coefficients, or Christoffel sym- 
bols, associated with the metric rjij (Chap. 4 in Ref. |28|). 
Introduce what is called the extrinsic curvature of the 
spatial surfaces as [s^ 



1 drji_ 



2a dC 

(see, e.g., Ref. or Sec. 21.5 in Ref. Then, 



a dt 



m7g 



2Ky + A', 



(77) 



(78) 



where Kj = rj^^Kij, in agreement with the standard rules 
of spatial index manipulation [Eqs. (I3ip ]. 

In Eq. ([78| , the first term on the right-hand side is due 
to the generally non-Euclidean form of r]ij , the second one 
is the gravity force, and the third one is due to the cur- 
vature of St considered as a subspace of spacetime (while 
r]ij , being the own metric of St , may or may not exhibit 
an intrinsic curvature) . The vector equation that we will 
now derive illustrates these forces in further detail. 



D. Vector equation 



It can be shown (Sec. 8.5 in Ref. [23|) that 



F 



kj 



dek 
dx^ 



(79) 



where are the basis three- vectors tangent to St, and 
are those of the dual three-vector basis. Then, 



1 df fyp" 



a dt 



+ 



k£ _ _ 

TO7 a 
_ 1 
a 



ek-^+P^{-v-V)ek 
dt 



djekp '') 
dt 



dek 



— e 



dc 



(80) 



The first term in the square brackets equals the fth pro- 
jection of dp/dt, whereas the second one can be repre- 
sented as follows. First, notice that the expression 



dek 
dC 



= e 



dek 

dc 



(81) 



coincides with the connection coefficient F^^ associated 
with the 4D metric g^i, (rather than r]ij), written in the 
basis B(. Similarly to Eq. ()76p . we can write then 

« _ fdgxc dgxk dgkc\ 
^ ''^ ~ ~ [~d^ ^ ~d^ W^)' 

where the underbars show that the expression is eval- 
uated in B(^. From Eq. (15^]) . both covariant and con- 
travariant metric elements with mixed (i.e., space-time) 
coefficients are zero. Then, since x^ = (, one gets 

j^e ^ ( ^9K 9gjk dgkQ \ 
2 \ dx^ dC dx3 ) 

_1 «%fe_l ^ary^ 

" 2 ^ dc ~ 2^^ dc ■ ^^^> 

Comparing this with Eq. (|77p. we obtain 
i dek 



dC 



-aKt 



(84) 



Hence, Eq. ([78]) rewrites as the following vector equation, 
independent of the spatial basis: 



1 dp 

a dt 



m7g -|-K-p + q(E + VxB). (85) 



[Note the unit coefficient in the term K • p, unlike in 
Eq. (j78l) .] For an alternative representation of Eq. ((85)) 
and comparison with similar representations found in lit- 
erature, see Appendix. 



IV. VLASOV EQUATION 

There are several ways to introduce the distribution 
function / describing the 3D motion in Sf that natu- 
rally extrapolate the one from the Minkowski spacetime. 
In this section, we will show that those definitions are 
equivalent; yet, the equation for / can take various forms. 
(Notice also that, when electromagnetic interactions are 
included, one may need to close the Vlasov theory with 
the Maxwell's equations. For those, see, e.g., Eqs. (3.4) 
in Ref. [Si] and also our Appendix for the notation.) 



A. Liouville theorem 

First, let us define / as the particle density in the 6D 
phase space, that is. 



/ = dN/d^, 



(86) 
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where dhf is the number of particles in the phase volume 
element dVl. In this case, the Vlasov equation can be 
derived immediately from the Liouville theorem (Sec. 3 in 
Ref. [S^l)- Namely, the latter says that dVt is conserved; 
then, since dAA is also constant, one gets 



dt 



0, 



(87) 



that is, / is conserved along the particle 6D trajectories. 

Yet these trajectories do not have to be expressed 
through canonical variables. For example, one may con- 
sider / as a function of and . Then, in agree- 
ment with Ref. [l5|, Eq. (|87l) rewrites like that for the 
Minkowski spacetime (see, e.g.. Chap. 8 of Ref. [HI), i.e., 

dl dx^dl d^dl^ 

dC dt dx' dt dp' ' ^ ' 

and one can use Eqs. ([57| and ([78]) for dx' /dt and dp^ /dt 
to close it. 



B. Divergence form 

To connect / with measurable quantities, dVt in 
Eq. (|86|) is derived as follows. From Sec. IIII Al it flows 
that the canonical variables are contravariant compo- 
nents and the canonical momenta are covariant com- 
ponents Pi. Since 



d{x\P,) 



= 1, 



(89) 



d{x^,Pj) 

one can also use kinetic momenta, though; then, 

dn^d^x*d^p^, (90) 

where d^x* = dx^ dx^ dx^ and cZ'^p* = dpidp2dp3 jssj. 
(From now on, asterisk denotes whether upper or lower 
indexes are assumed.) 

Hence, one needs to find how d'^x* and d^p* are con- 
nected with the invariant physical volumes in the coordi- 
nate and momentum spaces |14| . Since both of those are 
vector spaces with the metric rjij , one can write 



3 -* 



y/fjd-^x*, dVp = y/nd-^p 



Thus, 



dn = 



d{x\p,) 



dV^dVr, 



(91) 



(92) 



Since pi = rj^jP^ , the Jacobian here equals ry; therefore. 



dn = rfVxdVn 



(93) 



This means that /, originally defined as the phase-space 
density [Eq. ([55]) ] . is also the density in {x',p^) space: 



/ 



dAf 



dV^dVr, 



(94) 



in agreement with Ref. |14| . 

This result allows yet another representation of the 
Vlasov equation, which is derived as follows. First, 
rewrite Eq. as 



Vf 



dN 



d^x*d^p* 



(95) 



which means that rjf can be considered as the density /,j 
in {x^ ,p^) space, if the latter is assigned the Euclidean 
volume form (so the elementary volume is not ciVx^Vp 
but rather d^x* d^p*). Then, the particle current along 
i'-axis is = fr/dx' /dt, and the current along p*-axis 
is jp = frjdp' /dt, yielding that the particle conservation 
law reads as usual, 



dC, dx' dp' 







(96) 



(because the 6D space is considered Euclidean). Then, 
using that = rjf , one gets 



d f dx^ \ d (dp" 



(97) 



which we henceforth call the divergence form of the 
Vlasov equation. Equation (|97p. considered in combi- 
nation with Eq. for dx'/dt and Eq. ((TS]) for dp'/dt 
[and Eq. ([72]) for A'], represents the main result of this 
paper. In Sec.|Vl it will also be used to yield three- vector 
equations of coUisionless plasma hydrodynamics. 



C. Other representations 

The variables {x\p') are natural for describing dynam- 
ics on spatial hypersurfaces S*, because they allow for 
vector interpretation and are also self-contained (e.g., p'~' 
and po do not need to be considered). Yet, let us show 
how our formalism extrapolates to the original variables 
(x%p*), particularly, to compare with Ref. [l^ . 

In the form ((87|) . the Vlasov equation in variables 
{x\p') can be written immediately as 



df dx' df dp' df 
dt dt dx' dt dp' 



0, 



(98) 



in agreement with Ref. Correspondingly, the equa- 
tions for dx' /dt and dp' /dt can be obtained, e.g., from 
Eqs. dST]) and since 

dx^^dx^_nr d^^d^_ djP'pO) 
dt dt dt dt dt ^ ' 

[cf. Eqs. ([25]) and (g^])]. Here p° is given by 

dx° _ dt m7 

p = m — — = m — — mx = , (ILIU) 

dr dr a 
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T is the particle proper time, and [see Eq. (|52l) ] 

7 = \/l + riijp'^p' I'm? = -y/l + T^'-JpiPj/mP-, (101) 
the latter equality (to be used below) being due to 

p.,^p^. (102) 

Similarly, the divergence form is derived as follows. 
First, using Eq. (|102p . rewrite Eq. ([TO]) in the form 



d{x\pi) 



d{x3,Pj) 



(fx*d^p^ = (fx*d^p^ = rd^x*d^p*, 

(103) 



where r is the following Jacobian: 



dpi 



dpi 



(104) 



Combine Eq. (|103p with Eq. (|86l) . so one gets 

dM 
d^x*d'^p* 

Then, by analogy with Eq. ()97p . we immediately obtain 



(105) 



d{rf) 
dt 



d 

dx^ 



dx' \ d /dp' , 



(106) 



again in agreement with Ref . [iJI . 

Now let us show how the expression for r is derived 
(without introducing the "mass shell" used in Ref. p^). 
First, employ Eq. (|102l) . yielding 



opi dpi dpi 



(107) 



Then, using Eq. (|100p together with Eq. (|10ip . one gets 
dp^ /dpi = v' jo? , so Eq. (|107p rewrites as 



dpi 



n 



51 



-, (108) 



and therefore 1/r — w/rj, where w = det i. One of the 
ways to find the determinant of i, which is a 3 x 3 ma- 
trix, is through a brute-force calculation. A somewhat 
more elegant (and independent of the number of dimen- 
sions, albeit longer) way would be to see that w^i is a 
tensor of rank (1, 1), and thus w is independent of the 
(spatial) basis. Hence, one can consider the basis such 
that (3 points, say, along 2;-axis. Then, (3^ — 51(3^, in 
which case one immediately finds w = \ — VzP^ / ■ In 
the invariant form, this result is expressed through the 
scalar product of v and /3; therefore. 



1 - V /3/c 



(109) 



[When /3 = 0, one gets r — rj, and Eq. (|106p becomes 
equivalent to Eq. ([W)) . because (a;*,p*) are then the same 



as (x^p*) and t-axis coincides with ^-axis.] Using that 
Po — gOfiP^ , with g^jy expressed ^ in the basis Bq [rather 
than B(, as in Eq. (p9| ] and Eqs. ()100p and (|66|) . one can 
also rewrite Eq. (|109p as 



det g^j,^. 



Po 



Hence, our result agrees [3(| with that in Ref. \U. 



(110) 



V. HYDRODYNAMIC EQUATIONS 

Finally, let us consider moments of the Vlasov equation 
to obtain equations of coUisionless hydrodynamics. To 
preserve the vector form of the dynamic equations, we 
use (x*,p') variables. Correspondingly, the average of an 
arbitrary function x over the momentum distribution is 
defined as 



(111) 



where N stands for the particle density in the 3D space: 



N 



J fdVp^^J fd^p* 



(112) 



(not to be confused with the particle proper density [5s] ) . 
Integrating Eq. (|97|) over d'^p* yields 



d_ 

dc 



d 



(NVv) + (aNU^^) = 0, (113) 



where we introduced the flow velocity, as measured by 
FO, according to U = (V). Using the expression for the 
3D divergence operator written in the metric rjij (see, 
e.g.. Sec. 4.7 in Ref. 28]), one can rewrite Eq. (|113p as 



0. 



(114) 



This represents the continuity equation, which could also 
be obtained from the particle conservation in the 4D 
spacetime, by requiring that the four-divergence of the 
particle flow be zero |3| . 

As our next step, let us multiply Eq. (j97| by p> and 
then integrate over d^p* again. In this case, one obtains 

+ aNf{fi^V^ + [V • {ah)Y = aN{T + K • Vy . 

(115) 

Here we introduced the average momentum V ~ (p) , the 
average force on a particle 



T = m(7)g + K- P + g(E-fUxB), 



(116) 
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and the pressure tensor 

n^fc = f [pi - V1){V'' - U^) f dVp, 



(117) 



so one can interpret the fourth term in Eq. ()115p as jth 
component of its divergence (Sec. 4.7 in Ref. ^28,]): 



[V • {aU)Y 

With Eq. pi4p taken into account, Eq. pisp rewrites as 
1 dV^ 



a dC 



(K • vy + (u • v)v^ 
1 



= -^[v.(.n)F 



(119) 



Then, following the same argument as in Sec. IIIIDl one 
obtains the vector equation 



^__Lv.(an) 



+ m(7)g + K-7' + Q(E + U X B), (120) 

similar to that in the Minkowski metric. 

Higher moments of the Vlasov equation, which could 
yield a hydrodynamic closure like in Refs. [33-l4Q|, can be 
obtained analogously, and most easily for a nonrelativis- 
tic motion (7 ~ 1). Those are not discussed here, but 
are addressed separately in our Ref. [4l| . where we con- 
template the evolution of linear waves in a metric with 
nonzero K. Examples showing how the above equations 
can be used for applied calculations are also given in 
Ref. 



VI. CONCLUSIONS 

In this paper, the 3D dynamics of a charged particle 
in an arbitrary spacetime metric .g^,y, traditionally ad- 
dressed within differential geometry, is reformulated in 
terms of linear algebra (Sec. |Tl]) and Hamiltonian formal- 
ism (Sec. IIII[) . The modification of the Vlasov equation, 
in its standard form describing a charged particle distri- 
bution in the 6D phase space, is then derived explicitly, 
in two equivalent forms [Eqs. ([55)1 and (P7|]. The equa- 
tion accounts simultaneously for the Lorentz force and 
the effects of general relativity, with the latter appearing 
as the gravity force and an additional force due to the 
extrinsic curvature K of spatial hypersurfaces Et. For 
an arbitrary spatial metric, the equations of collision- 
less hydrodynamics are also obtained in the usual three- 
vector form [Eqs. (|114|) and (|120|) ]: for their applications, 
see our Ref. [41[. Another form of the Vlasov equation, 
which does not lead to vector equations but, on the other 
hand, allows for an arbitrary spacetime basis yBg, is also 
derived [Eqs. ((98|) and (|106|) ] within the new formalism 
and agrees with the results reported in Refs. [l^ [isj. 
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Appendix A: Another form of the particle motion 
equation 



Although Eq. ([85)) is sufficient for our purposes, let us 
explain how our notation relates to that from the widely 
cited Ref. Q (see also Refs. 0,11^). First of ah, the 
derivative on the left-hand side of Eq. ([55)) is the same 
as (the spatial part of) the four-vector 



h.to)=^-ih/^(h.p) 



dr-p 



dT 



(Al) 



where r is the particle proper time, dr = drpo/j; thus, 
Eq. (pp coincides with Eq. (3.8) in Ref. [3]- Following 
Ref. [3|, expand this derivative as 

= Z?,p + i (v • V)p = i?,p + (V • V)p. (A2) 

Here is called the Fermi- Walker derivative. 



1 dp 

a oQ 



(A3) 



which applies to both the components and the basis 
vectors that constitute the three- vector p = eip^. (The 
index r denotes differentiating with respect to Tpo rather 
than r.) One can further rewrite Z?rP as 



1 / dp^ dei 



Hence, its projection on the ith axis equals 

1 



(DrP) 



a dC 



(A4) 



(A5) 



where we used Eq. ([77]). Reverting to the vector form, 
one can write then 

DrP = VrP - K • p (A6) 

(cf. Eq. (2.16b) in Ref. |3!|), where, by definition. 



VrP = Vr{eip ) = T-r- 

a aQ 



(A7) 



That is, the new derivative Vr acts on the vector compo- 
nents but does not affect Gj. Notice now that, in the basis 
Bq, the partial derivatives of vector components with re- 
spect to C coincide with the Lie derivatives along C* = an 
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(see, e.g., Ref. @). Since this holds for aU components 
of I?rP: the latter can be given the following covariant 
definition: 

VrP = - {£a.pY, (A8) 

a 

that is, aVr can be understood as (the spatial part of) 
the Lie derivative along an (cf. Eq. (2.13) in Ref. [3j]). 

Similarly, one can introduce yet another Lie derivative, 
now along the original time axis Bq (rather than C*)- Us- 
ing Eq. (fTSj) . one obtains 

£t = £eo — £an + £ti- (^-9) 

Thus, for three-vectors, 

£tp = aVrP + £e,P (AlO) 

(cf. Eq. (2.16c) in Ref. [3]). Again use that, in the basis 
Bq, the Lie derivative along eo coincides with the partial 
derivative with respect to t: {£tpY — dp^/dt. Then, 

£,p^e,{£tpY ^e,^. (All) 

In other words, £t can be understood as the time deriva- 
tive, which differentiates the vector components but not 
the basis vectors in Bq. 

Combining the above formulas, one gets 

Dr = a-'{£t- £^)-K. (A12) 



Hence, the particle motion equation, Eq. (j85|) . rewrites 
equivalently as 



a M-ft + (v- V)]p = m7g + 2K-p + a ^£p,p 

+ q(E + V xB), (A13) 



where K can be alternatively put as 



K = -6--(6i/3)i, (A14) 



with (T being the traceless "shear tensor" , 6* = Tr K being 
the volume expansion rate, and I being the unit tensor. 
(Also, see Ref. [1] and Sec. 21.5 in Ref. [27] for the relation 
between K and n(a;'^).) This allows a direct comparison 
of our result with that given in Eq. (17) in the widely 
cited Ref. js^. Particularly, one can see that the latter 
must be corrected in the following aspects: (i) the time 
derivative that enters Eq. (jAlSp is, strictly speaking, the 
Lie derivative £t', (ii) the velocity v on the left-hand 
side is different from V on the right-hand side [Eq. ([56]) ] ; 
(iii) the Lorentz factor 7 is determined by V but not 
V [Eq. (155)) ]: (iv) the gravitational force mjg enters the 
right-hand side with the plus sign rather than the minus 
sign. 
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